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Abstract 



A wide range of new Q-conditional symmetries for reaction-diffusion systems with 
power diffusivities is constructed. The relevant non-Lie ansatze to reduce the reaction- 
PL,- diffusion systems to ODE systems and examples of exact solutions are obtained. Relation 

■ of the solutions obtained with the development of spatially inhomogeneous structures is 

discussed. 

PACS numbers: 02.20.-a, 02.30.Jr, 05.45.-a 

1 .Introduction. 



In 1952 A.C. Turing published the remarkable paper pQ, in which a revolutionary idea 
about mechanism of morphogenesis (the development of structures in an organism during the 
life) has been proposed. Roughly speaking, his idea says that the diffusion process can interact 
with the chemical reaction process in such a way that this can stimulate the development 



and growth of different forms and structures in an organism. Moreover this differentiation is 
often impossible without diffusion, i.e., only one may destabilize the spatially homogeneous 
structures. This effect is known as the Turing instability and mathematically leads to systems 
^ of reaction-diffusion equations with nonlinearities of the special form [21 [3]. Since 1952 the 

react ion- diffusion systems have been extensively studied by means of different mathematical 
methods, including group-theoretical methods. The main attention was paid to investigation 
of the two-component RD systems of the form 



U t = [D\U)U x } x + F(U,V) 
V t = [D 2 (V)V x ] x + G(U,V) 



2^u,i , n,TT^ (!) 



where U = U(t,x) and V = V(t,x) are two unknown functions representing the densities 
of cells and chemicals, F(U, V) and G(U, V) are two given functions describing interaction 
between them and environment, the functions -D 1 (f/) and D 2 (V) are the relevant diffusivities 
(usually they are assumed to be constants or power functions), and the subscripts t and x 
denote differentiation with respect to these variables. At the present time one can claim that 
all possible Lie symmetries of ([1]) with the constant diffusivities are completely described in 
[U El E] , while in [3 El [9] it has been done for the non-constant diffusivities. 

The problem of construction of conditional symmetries for (1) is still not solved even in the 
case of Q-conditional symmetries (non-classical symmetries). Moreover, to our best knowledge, 
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there are no any papers devoted to the search of conditional symmetries of the RD system (JTJ). 
Here, we present some recent results in this direction for the first time. It should be noted that 
there are many papers devoted to the construction of such symmetries for the scalar non-linear 
reaction-diffusion (RD) equations of the form g0\ QH CE21 13 HI [US [16] 

U t = [D(U)U x ] x + F(U). (2) 

and ([2]) with the convective term B(U)U X (here B(U), D(U) and F(U) are arbitrary smooth 
functions) [HI HSl EH] • 

It is well-known that conditional symmetries can be applied for finding exact solutions of 
the relevant equations, which are not obtainable by the classical Lie method. Moreover the 
solutions obtained in such a way may have a physical or biological interpretation (see, e.g., 
examples in [HJ [191 [SHI IZE] ) . In this paper we will demonstrate this using the Q-conditional 
symmetries of non- linear RD systems of the form ([T]). 

The paper is organized as follows. In the second section, we present new Q-conditional 
symmetries of the RD system 

U t = (U k U x ) x + F(U,V), 
V t = (V l V x ) x + G(U,V), 

which is the most important subcase of ([T]) from the applicability point of view [21 [22] • The 
determining equations for constructing Q-conditional symmetries of system ([3]) are derived and 
the main theorem presenting these symmetries in explicit form is proved. 

In the third section, the Q-conditional symmetries obtained are applied for reducing the 
corresponding RD systems to the systems of ordinary differential equations (ODE). A nonlinear 
RD system with the specified functions F(U, V) and G(U, V) is considered in details. We show 
that this system with the correctly-specified coefficients may lead to the Turing instability. 
Finally, we present some conclusions. 

2. Main result. 

It is well-known that the main difficulty arising in search of Q-conditional symmetries is 
to solve so called system of determining equations, which is nonlinear and overdetermined. 
Roughly speaking, one arrives on contradiction trying to construct Q-conditional symmetries 
of a given nonlinear PDE (system of PDEs) because there is the need to solve new nonlinear 
PDE system, which usually is much more cumbersome. This problem arises even in the case 
of linear PDE and it was the reason why G.Bluman and J.Cole in their pioneering work [23J 
were unable to describe all Q-conditional symmetries in explicit form even for the linear heat 
equation. In the recent paper [19], a short historical review is devoted to this problem and 
different kinds of non-Lie symmetry. 

Here, we construct the system of determining equations to find Q-conditional symmetry 
operators of the form 

Q = dt + £(t, x, U, V)d x + n\t, x, U, V)d v + V 2 (t, x, U, V)d v , (4) 

where £, rj 1 and rf are unknown functions, for the RD system Q and solve one under the rele- 
vant restrictions. Note that we search for purely conditional symmetry operators, which cannot 
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be reduced to Lie symmetry operators. We also remind the reader that each Q- conditional sym- 
metry operator is determined up to equivalent representations generated by multiplying (j3J on 
the arbitrary smooth function M(t,x, U, V). The problem of finding Q-conditional symmetry 
operators of the form 

q= e(t,x,^y)4 + r 7 1 (t,x,^y)a l/ + r / 2 (t,x,t/,y)ay, (5) 

is another challenge and will be treated elsewhere. 

First of the all, we remind the reader the notion of the conditional symmetry following the 
book [Till Section 5.7] (note one was formulated for the first time in |24|). 
Definition. A PDE of the form 

S (x 1 ,...,x r ,U,U,...u) = (6) 



(here U = U(xi, ...,x r ) and U is the totality of & th -order derivatives) is conditionally invariant 



under the operator 

Q = £\xi,..., x r , U)d Xl + ... + C( x i, %r, U)d Xr + rj(xi, x r , U)d v , (7) 

where rj and £ a with a = 1, . . . , r are smooth functions, if it is invariant (in Lie's sense) under 
this operator only together with an additional condition of the form 



S Q \x 1 ,...,x r ,U,y,...Uj=0, (8) 

that is, the overdetermined system of equations and (jSJ) is invariant under a Lie group 
generated by the operator Q. If the additional condition (JSJ) coincides with the equation Q(U) = 
0, i.e., 

^(Xi, X r , U)U X1 + ... + ^{Xi, X r , U)U Xr = 7)(xi, X r , U) 

then PDE ([6]) is Q- conditionally invariant under the operator (JTj). 

Obviously, this definition admits a direct generalization on systems of two PDEs. Of course, 
if the given system is not in involution then some difficulties can occur, however, system (pfl) 
has the involution form. Let us show how the determining equations to find Q-conditional 
symmetry operator (j3J) are obtained. 

It turns out that one can apply the local substitution 



u = U k+l , 
v = V l+l , I ^ -1. 



l+i i -/- 1 (9) 



to simplify the further computations. Of course, the case (k + 1)(Z + 1) = is special and needs 
of a separate investigation. Substitution © reduces system ([3]) and operator (j4]) to the form 



Uxx = U m U t + C 1 ^,^), 

Vxx = v n v t + C 2 (u,v), 



2 t „, „s (10) 



and 

Q = d t + ^(t,x,u,v)d x + r] 1 (t,x,u,v)d u + r] 2 (t,x,u,v)d v , (11) 

where m = ^ -1, n=-^j^-l, C 1 ^, i?) = -(k + 1)F (uw^, v^j , 

C\u, v) = -(p + 1)G (uw, , <9 n = sjj^-*^, d v = j^V- l d v . 
Let us apply now the second prolongation of the operator Q 



? 1 = Q+p 1 ti- t +Pti- t +p 1 ^+pi 



d 



~l~ < 7*<i; A,. . "I" "'(t fl„. ~t~ "t"Otta,„. ~^~^xx^. + C-r 



where coefficients p fc and cr fc with the relevant indices are calculated by the known formulas 
(see, e.g., P21I25]), to the each equations of (ITU1) . Since system (ITUj) should be considered as a 
manifold in the space of independent variables 

t, X, U, V, u t , v t , u x , V x , U xt , V xt , U t ti v u , u xx , v xx , 

we arrive at the invariance condition 

mn 1 u m ~ 1 u t + rfCl + rfC l v + p\u m = a\ x , 

nrj 2 v n-l Vt + V 2 C 2 + ^2 + p 2 v n = ^ l^J 

To obtain the determining equation for Q-conditional symmetry operator (1111) . one needs to 
take into account not only system fllOp (it will lead only to the determining equation for Lie 
symmetry operators !) but two additional conditions 

u t + £u x = r] 1 , v t + iv x = rf (13) 

generated by operator ffTTj) . Thus, inserting into ffT2l) the explicit expression for p k and a k , 
namely: 

P\ = Vt + Vl u t + vl V t ~ Uxfa + iuUt + CvVt), 

P 2 t = Vt + vl u t + Vv v t - v x (£ t + t,uU t + £>i), 

°lx = Vlx + 2 Vxu U x + 2 Vxv V x + vLi^x) 2 + vlviVxY + %vL U xVx + Vluxx + Vlvxx- 

Uxijixx ~\~ 2'Cxu'U'x ~\~ 2 £,xvA'x -\- ^uu(u x ) ~\~ ^vvi^x) ~\~ 2 ^uvU x V x -\- ^ U U XX -\- ^v^xx) (14) 
^l^xxx^x CuU x ~\~ ^v^x/i 

°L = Vlx + 2 vL u x + ZvLvx + vLM 2 + vlv^xf + 2 vL u xVx + vl^xx + r} 2 v v. 



XX 



Ux(£,xx 2 ^xu^x 2 Cxv^x £,uu(u x ) ~\- ^vviPx) 2 ^,uvU x V x ^u^ra ^v^xx) 
2 Vxx(^x ^u^x ~\~ ^v^x) 

and excluding four derivatives u t ,v t , u xx , v xx using (fTOj) and (fl3|) . one arrives at the cumbersome 
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expressions 



(15) 



u m (pi + rjfa 1 - fa + Tiltf - fa - u x + Uv 1 ~ fa) + Uv 2 ~ fa))) + 
+m V 1 u m - 1 (r ] 1 - fa + ^C* + ifCl = ^ + 2r)l u u x + 2r, l xv v x + vLfa? + vlfa 2 + 

^^Vuv^x^x U x (Jl xx ~\- 2^ XU U X -\~ 2^ XV V X -\~ £,uu{u x ) ~\~ C,vv{u x ) -\~ 2^ uv U x V x )-\- 

+ ((V 1 - fa> m + C 1 ) (rji - 2& - 3£ u u x - 2^) + (fa 2 - fav n + C 2 ) (^ - £ v u x ) , 

v n (vt + vKv 1 - fa) + v 2 v (v 2 - fa) - v x [it + Uv 1 - fa) + Uv 2 - fa))) + 

+nr ] 2 v n '\r ] 2 - fa + tfC 2 + V 2 C 2 = rf xx + 2rg u u x + 2rf xv v x + vLfa) 2 + vlfa 2 + 

^^Vuv^x^z V x ^ xx -\- 2^ XU U X -\- 2^ XV V X -\- ^uu{u x ) -\- ^vv{u x ) ~\~ 2^ uv U x V x )-\- 

+ (V - fafa + C 2 ) { V l - 2£ x - 3£ v v x - 2£ A )+ (fa 1 - iu x )u m + C 1 ) fa* - £ u v x ) . 

Finally, we take into account that the unknown functions r] 1 , rf and £ don't 
depend on the derivatives u x and v x and therefore we split two expressions 
arising in ([15]) on fa) 3 , u x (v x ) 2 , v x (u x ) 2 , u x v x , {u x ) 2 , (v x ) 2 , v x , u x and 
fa 3 , v x (u x ) 2 , u x (v x ) 2 , u x v x , fa) 2 , fa 2 , v x , u x , respectively. Thus, we arrive at the 
nonlinear system of determining equations 

1) £,vv £,uv 0, 

2) vL = o, 

3) vL = o, 

4) 2^ u u m + rjl u - 2U = 0, 

5) 2^ + ^-26^ = 0, 

6) Zfa m + v n ) + 2 V l-2Z xv = 0, 

7) ^fa m + v n ) + 2i 1 2 uv -2U = ^ 

8) fa(u m - v n ) + 2r,l v - 2^fa - 2fau m = 



9) fa(v n - u m ) + 2rg u - 2i u C 2 - = 0, (16) 

10) - mfau m - 1 + {2 fa -i t - fa ~ 2fa)u m + 
+fa 2 v n + 3 fa 1 + e,C 2 - 2r\\ u + i xx = 0, 

11) - nfav n ~ l + {2ifa - 6 - fa 1 - 2& x )v n + 
+fau m + 3 fa 2 + &C 1 - 2rg v + U = 0, 

12) mfafu™- 1 + (r)l + rfril + 2fa)u m - r] 2 r]fa+ 
HVC* + rfCl - viC 1 + 2Z tX C 1 - V IC 2 - r£. = 0, 

13) n(rf) 2 v n ~ l + {rj 2 + rfrg + 2i x rf)v n - r l 1 rf u u m + 
WC 2 U + v 2 C 2 - vlC 1 + 2i x C 2 - V 2 C 2 - rf xx = 0. 



to find the coefficients of the operator (TTTj) and the functions C 1 ,C 2 '. Equations 1) from this 
system are easily integrated and lead to 

£ = a(t, x)u + b(t, x)v + c(t, x), (17) 

where a, b and c are arbitrary (at the moment) smooth functions. Substituting (JTTJ) into equa- 
tions 6) and 7) of ffT6l) and taking into account the second and third equations of (|16p . one arrives 
at the requirement a = b = 0. Thus, equations 2)-7) of system (fT6l) can be straightforward 
integrated and their general solution takes the form 

f = c(t,x), 

t? 1 = q l (t)v + r l (t,x)u + p 1 (t,x), (18) 
rf = q 2 (t)u + r 2 (t, x)v + p 2 (t, x), 

where the functions in right-hand-side are arbitrary those of their arguments. 

The remaining equations 8)-13) of system (fToT) involve the unknown functions C 1 and C 2 
and are called the classification equations. To solve them one should consider three different 
cases depending on the functions q x {t), q 2 (t) and £(t, x) arising in (ITS]) : 

a) q\t) = q 2 (t) = 0, £(t,z)^0; 

b) q\t) = q 2 (t) = 0, f(t,a;) = 0; 

c) gH^ + ^WVO, e(^) = 0. 

Note that the fourth possible case g 1 ^) 2 + g 2 (£) 2 ^ 0, £(t,a;) 7^ arises only under the 
restriction m = n = 0, which follows from 8)-9) of system ffTB]) . It means the case of the 
RD system ([3]) with constant diffusivities should be separately considered so that we assume 
m 2 + n 2 ^ below. 

Solving equations 8)-13) of system (fT6|) in the case (a) is rather simple but very cumbersome. 
After the relevant computations ( the program package MATHEMATICA 5.0 was also used) 
it has been established that all the operators obtained are nothing else but the Lie symmetry 
operators (up to the relevant multiplier M(t, x, U, V)) found earlier in [S]. In contrary to (a), 
the case (c) is the most difficult and at the present time we were able to solve equations 8)-13) 
of system f|T6l) only in particular subcases. 

Now we present the theorem giving a complete description of Q-conditional symmetry op- 
erators of the form (J3J) in the case (b). 

Theorem 1 The RD system ^ with k 2 + I 2 ^ and [k + + 1) 7^ is Q- conditional 
invariant under the operator with £ = and rjy = rjfj = if and only if it and the relevant 
operator (up to equivalent representations generated by multiplying on the arbitrary smooth 
function M(t, x, U, V) ) have the forms listed in the table 1 (here f and g are arbitrary smooth 
functions of the relevant argument, while Xj,j = 1,2,3,4 are arbitrary constants). 



Table 1. Q-conditional symmetries of the RD system ((3D 



no 


RD systems of the form ([]} 


Q-conditional operators 


Restrictions 


1. 


U t = (U k U x ) x + f{U k+1 ), 
V t = (V-iv x ) x - 2\V? + g(U k+1 ) 


d t + 2p{x)V?d v 


Pxx = (P) + Xp, 

P^O 


2. 


U t = (U k U x ) x + MU- k + f(U k+1 - aV l+1 ), 
V t = (V l V x ) x + X 2 V- 1 + g{U k+1 - aV l+1 ) 


d t + XiU- k du+ 

+\ 2 v- l d v 


_ Ai(fc+i) X / n 
" - A 2 (i+1) ' 2 ' " 

A? + Z V 0. 


3. 


U t = (U~^U X ) X - 2XU^ + /(C/3 - V*), 
V t = (V~iV x ) x - 2XVi + g{ui ~ V*) 


d t + 

+2p(x)(Uid u + Vid v ) 


Pzz = (J?) 2 + Ap, 

P^O 


4. 


U t = (U k U x ) x + \ 1 U- k + /(w), 
V t = (W x ) x + (V l+1 - A 3 )(.gM + \ 2 V- 1 ) 


d t + XiU- k du+ 
+\ 2 {V - \ 3 V- l )d v 


w= ^^Un, A 2 ^0, 

(yi + i_A 3 ) ^2<' + i) 

eii/ier A 2 + A§ 7^ or 
A| + k 2 ^ orA 2 + IV 


5. 


u t = (u k u x ) x + (u k+1 - Ai)(/(w) + \ 2 u- k ) 

V = (V l V x ) x + (V l+1 - \ 3 )(g{u) + X 4 V~ l ) 


dt + X 2 (U-X 1 U- k )d u + 
+X 4 (V - X 3 V- l )d v 


w = x/(fc+i) i A2A4 7^ 0, 

eii/ier A 2 + A 2 7^ 0, 
or X 2 , + k 2 ^ or X\ + I 2 ^ 



Sketch of the proof of Theorem. To prove the theorem one needs to construct the 
general solution of subsystem 8)-13) of system (JIB]) under the assumption £ = and r?y = 

= 0. Obviously, equations 8) and 9) are automatically satisfied, while those 10) and 11) are 
reduced to the form rv\ v = and rj 2 v = 0, respectively, i.e.: 



r^r 1 ^), r 2 = r 2 (t). 
So, the last equations 12) and 13) take the form 
( r i M + pi)c£ + ( r 2 t; + p 2 )^ 1 - 

+ (rj + m^ 1 ) 2 )^ 1 + (p] + 2mr 1 p 1 )u m + m^fvJ 11 ' 1 - p^, = 0, 
(r 1 ^ + p v )C 2 u + (r 2 w + p 2 )Cl - r 2 C 2 + 

+ (r 2 + n(r 2 ) 2 )v n+1 + (p 2 + 2nr 2 p 2 )v n + n(p 2 ) 2 7j n_1 - p 2 ^ = 0. 

System (|T9|) consists of two independent first-order linear PDE with respect to the unknown 
functions C l (u, v) and C 2 (u, v) therefore its general solution can be straightforward constructed, 
however we should remember that the coefficient in (I19p are functions on t and x. To construct 
all possible solutions of ffl~9]) one needs to consider six different cases (up to renaming u — > v 
and v —> u): 

1) r 1 = p 1 = r 2 = p 2 = 0, 

2) r 1 = p 1 = 0, r 2 ^ 0, 

3) r 1 = p 1 = r 2 = 0, p 2 ^ 0, 

4) r 1 = 0, pV 0, r 2 = 0, p 2 ^ 0, 

5) r 1 = 0, pV 0, r 2 ^ 0, 

6) r 1 ^ 0, r 2 ^ 0. 
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In the case 1) operator (j3J) immediately takes the form Q = d t , which is, of course, the Lie 
symmetry operator. The similar situation occurs in the case 2) since all the operators obtained 
are equivalent to the relevant Lie symmetry operators listed in [Sj. The most interesting cases 
are those 3)-6). 

Consider case 3) in details. In this case system ffT9l) takes the form 

P C v = 0, , , 

p 2 C 2 + p 2 t v n + n{p 2 fv n - 1 - p 2 xx = [Z()) 



and its formal integration leads to the solution 
C 1 = f(u) 

C 2 = l(^~ - npV -1 ) dv + g(u), 



(21) 



where / and g are arbitrary smooth functions. Since the function C 2 doesn't depend on on t 
and x, three subcases should be separately examined: n = 0, n = 1 and n ^ 0; 1. The first 
subcase immediately gives C 2 = + g(u), so that 

Pxx ~ Pt _ \ 
2 _ A ' 

where A is an arbitrary constant. So, the system 

u xx = u m u t + f(u), , 22 s 
v xx = v t + Xv + g(u) 

admits the Q-conditional symmetry operator 

Q = d t +p 2 (t,x)d v , (23) 

where p 2 (t,x) is the general solution of the linear PDE p 2 t = p xx — Xp 2 . However, if one now 
applies substitution (J9j) to f[22l and (|23|) . then the RD system and the operator listed in [8] 
(see case 5 in table 1) are obtained. So, subcase n = doesn't lead to any Q-conditional 
symmetries. 

In the subcase n = 1 the general solution of (120]) takes the form 
C 1 = /(«), C 2 = Xv +g(u), 

2 

where A = ^S- — p 2 . So, the system 

u xx = u rn u t + f(u), , 24 , 
Vxx = vv t + Xv + g{u) " 

admits the Q-conditional symmetry operator 

Q = d t +p 2 (x)d V} (25) 
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where the function p 2 (x) is the general solution of the nonlinear ODE 

PL = (P 2 ) 2 + Ap 2 . (26) 

Applying now substitution to (|24l) - (|25j) and introducing the relevant notations, one arrives 
at the system and the Q-conditional operator listed in case 1 of table 1. 

Considering the subcase n ^ 0; 1, we immediately obtain p 2 = A = const (see ( 121]) ) and this 
lead to the system 

u xx = u m u t + f(u), 
v xx = v n v t - Xv n + g{u) 

and the operator 

Q = d t + \d v . (28) 

Operator f[2"gj) is reduced to the form Q = d t + ■^■XV~ l dy with A ^ by using substitution (J9j). 
On the other hand, system ( 1271) and operator ( 1251) correspond to a particular case at Ai = a = 
of those listed in case 2 of table 1. Thus, case 3) is completely investigated. 
The case 4) can be examined in a quite similar way and the system 

u xx = u m u t — aXu m + f{u — av), 

v xx = v n v t - Xv n + g(u - av), (29) 
Q = d t + X(ad u + d v ), 

and the operator 

Q = d t + X(ad u + d v ) (30) 

are obtained, where a ^ is an arbitrary constant. It is easily seen that systems and operators 
(l2~7j) - (!30l can be united, i.e., the restriction a ^ is not essential. Applying now substitution 
([9]) to (!2~9l) - (!30l) and introducing the relevant notations, one arrives at the system and the Q- 
conditional operator listed in case 2 of table 1. It turns out that the power m = n = 1 leads to 
an additional symmetry in this case. In fact, the system 

u xx = uu t + Xu + f(u - v), , . 

v xx = vv t + Xv + g(u - v) 

is conditionally invariant with respect to the operator 

Q = d t +p 2 (x)(d u + d v ) } (32) 

where p 2 (x) is the general solution of the nonlinear ODE (|26|) . Formulas (!31~!) - (!32l together 
with the substitution ([9]) generate the RD system and the operator listed in case 3 of table 1. 
Those listed in cases 4 and 5 of the table can be similarly obtained by examination of the cases 
5) r 1 = 0, 0, rV and 6) r 1 ^0, r 2 ^ 0. 

Finally, we note that all operators arising in table 1 are not Lie symmetry operators because 
any Lie symmetry operator of the RD system ([3]) must be linear on U and V [8j. 

The sketch of the proof is now completed. 

Remark 1. Restrictions on the coefficients A^ arising in the last column of table 1 guarantee 
that the relevant operators do not coincide with Lie symmetries. Of course, those operators is 
still Q-conditional symmetry operator if some lambda-s vanish and then they are equivalent to 
the relevant Lie symmetry operators obtained in [S]. 
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3. Exact solutions and their application. 



In this section, we apply the operators of conditional symmetry for constructing exact 
solutions of the relevant RD systems, demonstrate their application for solving boundary value 
problems and present some biological interpretation. 

It is well-known that any Q-conditional symmetry operator of the given two-dimensional 
PDE (system of PDEs) guarantees its reduction to an ODE (system of ODEs) by construction 
of the corresponding local substitution, which is called ansatz. Usually the ODE (system of 
ODEs) obtained can be integrated (at least partly) and exact solutions of the initial PDE 
(system of PDEs) are constructed. Nevertheless, this procedure contains new difficulties, one 
can successfully realize one for the operators arising in table 1. We present the final result in 
the form of table 2. 

Table 2. Ansatze and reduced systems of ODEs for RD systems of the form ((3]) 



no 


Ansatze 


Systems of ODEs 


1. 


U = (<p(x))T^ 
V = (p 2 t + 4>(x)) 2 


fxx = ~(k + l)f(<p) 

^xx = -\g{v) + (A + p)i/j 


2. 


U = (X 1 (k + l)t + <f(x))T^ 

V= (X 2 (l + l)t + lfj(x))TTT 


Vxx = ~{k + l)f((p - aip) 

i>xx = -(I + i)g{<p - oiiji) 


3. 


U = (p 2 t + ip(x)) 2 

v = ( P H + ip{x)f 


tPxx = (p+X)lf-p(lf-1p) 

ipxx = (p + X)tjj - \g(ip - ip) 


4. 


U= (X^k + l)t + <p(x))T& 
V = (exp(A 2 (/ + l)t + ip(x)) + X 3 )ttt 


<p xx = (k + 1)/ (exp(<P a^+I)^)) 
i>xx + € = -{I + 1)9 (exp(v^ - j^Sj^)) 


5. 


U = (cp(x) exp(X 2 (k + l)t) + Ai)^t 
V = (ip(x) exp(A 4 (Z + l)t) + A 3 )^ 


tp xx = -(k + l)tpf UiP Mm) J 

/ A 2 (fc+1)\ 

ipxx = -(l + l)ipg Upip A ^ +1 > J 



Remark 2. The numbers 1,. . . ,5 in table 2 correspond to the RD systems arising in table 
1 with the same numbers. 

Now one sees that the reduced systems of ODEs are nonlinear and it is quite implausible 
that those are integrable for arbitrary smooth functions / and g. However, these systems can 
be integrated if the functions / and g are correctly specified. For example, the ODE system 
arising in case 5 takes the form 

<Pxx = aiV + M, ( 33 n 



in 



if one sets g = -^(a 2 J + / = -^("1 + A J), r = -A a (A; + 1) = -A 4 (Z + 1), where 
ctj, fy, j — 1,2 are arbitrary constants. Nevertheless the initial RD system 

U t = (U k U x ) x -^{rU + ai (U k+1 - Ai) + - A 3 ) - XirU-% 

V t = (V l V x ) x - j^(rV + a 2 (U k+1 - X,) + p 2 (V l+1 - A 3 ) - X 3 rV~ l ) 

is still non-linear, the reducing system (155]) is linear and its general solution can be constructed 
in explicit form. Assuming (3\ ^ (otherwise should be a 2 ^ and one will start from the 
second equation of (1551) ) the function ^ can be expressed from the first equation of (1551) so that 
the second equation takes the form 

fxxaix ~ («i + P2)¥xx + (ai/?2 - a 2 Pi)<p = 0. (35) 

The general solution of the fourth order ODE (155]) essentially depends on the roots of the 
algebraic equation 

s 4 - («! + (3 2 )s 2 + (ai/% - a 2 Pi) = 0. (36) 

Generally speaking, 9 different forms of the function (p can be obtained. To avoid cumbersome 
computations, we consider only the case when (156]) possesses four different complex roots with 
the zero real parts, i.e.: 



'-(ai+ft + i/Oi- 


-ft>) 2 +4a 2 /3i) 


2 


-(ai+fe-yfar 


-/3 2 ) 2 +4a 2 /3i) 



si >2 = ±i 1 

It occurs under the following restrictions on the coefficient of (1551) : 



(37) 



(«! -/9 2 ) 2 > -4a 2 /?i, ai + & < - \/(«i - Z^) 2 + 4a 2 0i- (3* 
Hence, we obtain the general solution of (I55j) : 



9? = Ai cos(\si\x) + A 2 sin(|si|x) + A 3 cos(|s3|x) + sin(|s3|x), 
if; = — ^ ^(ai + |si| 2 )(v4i cos(|si|x) + A 2 sin(|si|x)) + 



(39) 



+{oci + |s3| 2 )(A 3 cos(|s3|2;) + y4 4 sin(|s 3 |x)) 
where A 3 j = 1,2,3,4 are arbitrary constants. Substituting (159"|) into ansatz 



U = (<p(x) exp(X 2 (k + l)t) + Ai)*+i 
V = (^(a;)exp(A4(/ + l)t) + A3)^T, 



(40) 



arising in case 5 of table 2, we arrive at the fourth-parametrical family of exact solutions of 
the nonlinear RD system (1541) with the coefficients satisfying restrictions (158]) . It should be 
noted that all solutions of the form (1591 - (140]) are periodic on the space variable x. According 



1 1 



to J.D. Murray [2], the periodic (in space) solutions can mathematically express the Turing 
instability, which leads to the spatially inhomogeneous structures and forms. Moreover, we can 
show that solutions of the form (p39l) - (T4T3l) with the correctly specified constants Aj,j = 1 . . .4 
satisfy the zero Neumann conditions (zero flux boundary conditions), which naturally arise in 
mathematical models with the Turing instability. In fact, consider the space interval [0, a], a > 
where the cells with the density U and chemicals with the density V are interacting. Assuming 
that their interaction is described by the RD system (1511) . the zero flux boundary conditions 
on [0, a] lead to the requirements 

U x \ x =o = 0, V^l^o = 0, , . - 

U I =0 V I =0 ^ ' 

LJ x\x=a w ; y x\x=a 

Using ansatz (HOI) , one obtains formulas 

k 

U x = m ip x (ipexp(-rt) + Ai)~*+i, 



V x = ITI^O exp(-rt) + A 3 ) i+i 



Hence, taking into account solution (I39p and conditions (j41j) with x = 0, we arrive at two linear 
algebraic equations 

<fx\x=o = A 2 \si\ + A 4 \s 3 \ = 0, 

ip x \x=o = A 2 \si\(a 1 + |si| 2 ) + A 4 |s 3 |(ai + |s 3 | 2 ) = 

to find Ai and A4. Analyzing these equations one easily finds that the solution A 2 = A 4 = 
only is possible. Finally, using conditions (l4~Tj) with x = a, we obtain 

sin(|si|a) + A 3 \s 3 \ sin(|s 3 |a) = 0, . . 

^4i|si|(ai + |si| 2 ) sin(|si|a) + A 3 |s 3 |(ai + |s 3 | 2 ) sin(|s 3 |a) = 

to find A\ and A3. Analyzing (j4"51) . we deduce that the space interval [0,a], a > cannot be 
arbitrary otherwise tp = ip = 0. There are three countable sets of values for a: (i) a — j^r, 

(ii) a — and (in) a = 7 L^EM2= : (ji, J2) C N 2 , for which condition (14"5|) is satisfied by the 
functions 

ip = At cos(|si|x), ip — ~ -?rAi(a\ + |si| 2 ) cos(|si|x), 

Pi 

<p = A 3 cos(|s 3 |x), ip = -—A 3 (ai + \s 3 \ 2 ) cos(|s 3 |x), 

Pi 



and 



tp = Ai cos(|si|x) + A 3 cos(^-\si\x) , 



ip = -^(A 1 (a 1 + |si| 2 ) cos(|si|x) + A 3 (a 1 + j||si| 2 ) cos(g|si|x)^ . 
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respectively. Moreover, the additional restriction «2 = (ai ~' 2 ^^jf^^ ai31 ^ is obtained in the 
case (Hi), which allows to simplify the expression for \sx\ and | S3 1 : 



si = -^ . 9 s — , S3 = — si • (44) 

Thus, we can now formulate the theorem. 

Theorem 2 The nonlinear RD system ft34\ ) u>zi/i the coefficients satisfying restrictions 
possesses the periodic (in space) solutions 



1 

fc+i 



[i) U =( A% cos(|si|x) exp(— rt) + Ai 

1 ^ = I - (Ax(ax + \sx\ 2 ) cos(\sx\x)) e ^f) + x 3 



1 

i+i 



(45) 



1 

fc+i 



satisfying the zero boundary conditions on the interval 
(ii) U = ( A 3 cos(|s 3 |a;) exp(— rt) + Ai 

\ ' = I - (A 3 (ax + |s 3 | 2 ) COS (\s 3 \x)) c -^^ + A 3 



0,pt 



, h e N; 



1 



satisfying these conditions on the interval 



' l«3| 



(w) Z7 =( (Ax cos(|si|x) + A 3 cos(4^|si|x)) exp(— rt) + Ai 



1 

fc+i 



(46) 



V = ( - ( Ax(ax + \sx\ 2 ) cos(\sx\x) + A 3 (ax + ||si| 2 ) cos(g|si|x) ] , 



j£ l. 1 2^ „„ D cii| D j^ ^ cx p(- rf ) _|_ \ 3 N 



j_ (47) 



satisfying these conditions on the interval 



0. 



, 0i,J 2 )cN 2 



cas 



v / -( Q i+/ 3 2) 

In t/ie cases (i)and (ii), the values of \sx\ and \s 3 \ are determined by formulae pT7| ) ; m i/ie 

(aij2-feJi)(feJ2- Q iJ?) 

(ii 2 +i 2 2 ) 2 A ' 



e fraj t/ie uaZue of \s\\ is given in fiJ4\ ) and cii 



Remark 3. Since system (151)1 and boundary conditions (|4"T1) are invariant under the space 
translation transformations where x$ is an arbitrary parameter, one can generalize 

this theorem on the intervals [— a, 0] and [—a, a]. 

One sees that each solution of the forms (j4"5]) - (j4"T|) tends to the steady state point (U, V) = 
((Ai) ~ , (X^ 1 ^) if the time t — > +00 and r > 0. This point can be stable or unstable depending 
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Figure 1: Exact solution (|45]) with r = 2, A x = 0.95, s x = \Jl - V2, A x = 1, A 3 = 2 




Figure 2: Exact solution gS} with r = -2, Ay = 0.1, Si = a/2 - \/2, A x = 1, A 3 = 2 



on the coefficients of system An example of solution ( H5l) of the RD system 

K = (wg. - 1^ + f 72 + v 2 + ^ (48) 

with with r = 2 is presented on Fig.l. This solution arrives at the steady state point (U, V) = 
(1, \/2) very fast. 

Solutions (H5|) - (|47|) with r < are growing to infinity with time providing k + 1 > and 
I + 1 > 0. On the other hand, the spatially inhomogeneous structures should be formed for a 
finite time (the development of those structures in an organism should be finished during the 
organism life). Taking this into account, one notes that solutions fH5l) - fH7j) with r < may 
describe the development of spatially inhomogeneous structures for the finite time t max . The 
relevant example of solution fH5l) of the RD system fHHl) with r = —2 is presented on Fig.2. 

One sees that both components, U and V have practically homogeneous distributions in 
space if t « t max . However, those are essentially inhomogeneous if t » 0. Moreover each 
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of the components U and V dominates in the different regions. So, we may say that this 
solution forms a striped structure. This kind of spatially inhomogeneous structures is typical 
for mammalian coat (tigers, zebras). 

4. Conclusions. 

In this paper, theorem 1 giving a complete description of Q-conditional symmetries of the 
nonlinear RD system ([3D under some restrictions is proved. We established that there is a 
wide range of the nonlinear RD systems of five types, which are conditionally invariant with 
respect to the operator PJ with the coefficients £ = and Vv = Vu = 0- Moreover those 
systems don't coincide with the RD systems admitting a non-trivial Lie symmetry, which has 
been completely described in [3 [8]. To our best knowledge, there are no any papers devoted to 
constructions of conditional symmetries for systems of evolution equations, excepting the recent 
paper [26]. Note that only an example of the nonlinear diffusion-convection system was found 
in [26] , while here we present five classes of the nonlinear RD systems and each of them involves 
two arbitrary functions. The work is in progress to solve the problem without any restrictions 
on the coefficients of operator ([3]). 

The Q-conditional symmetry operators have been applied for obtaining exact solutions of 
the relevant RD systems. In fact, we reduced the problem of solving PDE systems to one for 
the relevant ODE systems. These ODE systems were further solved therefore several exact 
solutions in explicit form have been constructed for the RD systems (}34"j) . Properties of the 
solutions obtained were investigated and a possible biological interpretation was suggested. 

Acknowledgments 

The authors are grateful to the unknown board member and the unknown referee for the 
useful comments. 

References 

[1] Turing A M 1952 The chemical basis of morphogenesis Phil. Trans. Roy. Soc. Lond. 237 
37-72 

[2] Murray J D 1989 Mathematical Biology (Berlin: Springer) 

[3] Britton N F 2003 Essential Mathematical Biology (Berlin: Springer) 

[4] Cherniha R, King J R 2000 Lie Symmetries of Nonlinear Multidimensional Reaction- 
Diffusion Systems:I J. Phys. A:Math. Gen. 33 267-282, 7839-41 

[5] Cherniha R, King J R 2003 Lie Symmetries of Nonlinear Multidimensional Reaction- 
Diffusion Systems:II J. Phys. A: Math. Gen. 36 405-425 

[6] Nikitin A G 2005 Group classification of systems of non-linear react ion- diffusion equations 
Ukrainian Math. Bull. 2 153-204 

[7] Knyazeva I V, Popov M D 1994 A system of two diffusion equations, in: Lie Group Analysis 
of Differential Equations, Editor N.H. Ibragimov. Boca Raton: CRC Press 1 171-176 

15 



[8] Cherniha R, King J R 2005 Nonlinear Reaction-Diffusion Systems with Variable Diffusiv- 
ities: Lie Symmetries, Ansatze and Exact Solutions. J. Math. Anal. Appl. 308 f 1-35 

[9] Cherniha R and King J R 2006 Lie Symmetries and Conservation Laws of Nonlinear 
Multidimensional Reaction-Diffusion Systems with Variable Diffusivities, IMA J. Appl. 
Math. 71 391-408 

[10] Fushchych W I, Shtelen W M, Serov M I 1993 Symmetry analysis and exact solutions of 
equations of nonlinear mathematical physics (Dordrecht: Kluwer) 

[11] Serov M I 1990 Conditional invariance and exact solutions of non-linear heat equation 
Ukrainian Math. J. 42 1370-76 

[12] Nucci M C 1992 Symmetries of linear, C-integrable, S'-integrable and nonintegrable equa- 
tions and dynamical systems Nonlinear evolution equations and dynamical systems (River 
Edge: World Sci. Publ. (USA)) 374-381 

[13] Clarkson P A and Mansfield E L 1993 Symmetry reductions and exact solutions of a class 
of nonlinear heat equations Physica D 70 250-288 

[14] Arrigo D J, Broadbridge P and Hill J M 1994 Nonclassical symmetry reductions of the 
linear diffusion equation with a nonlinear source. IMA J. Appl. Math. 52 1-24. 

[15] Arrigo D J, Hill J M 1995 Nonclassical symmetries for nonlinear diffusion and absorption 
Stud. Appl. Math. 94 21-39 

[16] Pucci E, Saccomandi G 2000 Evolution equations, invariant surface conditions and func- 
tional separation of variables Phys. D 139 28-47 

[17] Cherniha R, Serov M 1998 Symmetries, Ansatze and Exact Solutions of Nonlinear Second- 
order Evolution Equations with Convection Term Euro. J. Appl. Math. 9 527-542 

[18] Cherniha R 2007 New Q-conditional Symmetries and Exact Solutions of Some Reaction- 
Diffusion-Convection Equations Arising in Mathematical Biology J. Math. Anal. Appl. 326 
783-799 

[19] Cherniha R, Pliukhin O 2007 New conditional symmetries and exact solutions of nonlinear 
reaction-diffusion-convection equations. J. Phys. A 40, 10049-70. 

[20] Dixon J M, Tuszynski J A, Clarkson P A 1997 From Nonlinearity to Coherence (Oxford: 
Clarendon Press) 

[21] Bradshaw-Hajek B H, Broadbridge P 2004 A robust cubic reaction-diffusion system for 
gene propagation. Math. Comput. Modelling 39, no. 9-10 1151-63. 

[22] Okubo A. and Levin S.A. 2001. Diffusion and Ecological Problems. Modern Perspectives. 
2-nd erf. (Berlin: Springer) 



1fi 



[23] Bluman G W and Cole J D 1969 The general similarity solution of the heat equation 
J. Math. Mech. 18 1025-42 

[24] Fushchych W I, Serov M I, Chopyk W 1988 Conditional invariance and nonlinear heat 
equations Dopovidi Akad.Nauk Ukrainy Ser. A (Proc. Ukrainian Acad. Sci. Ser. A) 9 
17-21 (in Ukrainian) 

[25] Olver P. 1986 Applications of Lie Groups to Differential Equations ( Berlin: Springer) 

[26] Cherniha R, Serov M 2003 Nonlinear Systems of the Burgers-type Equations: Lie and Q- 
conditional Symmetries, Ansatze and Solutions J.Math.Anal.Appl. 282 305-328. 



17 



